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Summary. We present a study of a peculiar Supernova Remnant MF16, associated 
with the Ultraluminous X-ray Source (ULX) NGC6946 ULX-1. Observations were 
taken with the MultiPupil Fiber Spectrograph (MPFS) with the 6-m telescope on 
January 2005. The nebula is found to be highly asymmetric, one of the parts being 
much denser and colder. Two-component structure of the emission lines and radial 
velocity gradient in some of them argue for a non-spherical nebula, expanding with 
a velocity of about lOOfcms"^. Neither shock models nor the X-ray emission can 
adequately explain the actual emission line spectrum of MF16, so we suggest an 
additional ultraviolet source with a luminosity of about lO^^erg s~^. We confirm 
coincidence of the ULX with the central star, and identify radio emission observed 
by VLA with the most dense part of the nebula. 



Introduction 

The attention to MF16 was first drawn by Blair & Fesen [1], who identified 
the object as a Supernova Remnant (SNR), according to the emission-line 
spectrum with bright coUisionally-excited lines. It was long considered as a 
peculiar and unusually luminous SNR, because of its huge optical emission- 
line {Lho ~ 1-9 X lO^^ergs^^, according to [1], for the tangential size 20 x 
34pc) and X-ray {Lx ~ 2.5 x lO'^^erg in the 0.5 — 8keV range, according 
to [2]) luminosities. 

However, it was shown by Roberts & Colbert [2] that the spectral, spatial 
and timing properties of the X-ray source do not agree with the suggestion 
of a bright SNR, but rather suppose a point source with a typical "ULX- 
like" X-ray spectrum: cool Multicolor Disk (MCD) and a Power Law (PL) 
component. So, apart from the physical nature of the object, MF16 should be 
considered a ULX nebula, one of a new class of objects, described by PakuU 
& Mirioni, [3]. For introduction into ULX-counterparts see [4]. 



Observations and Data Reduction 

We analyse spectra obtained with the SAO 6m telescope using Multipupil 
Fiber Spectrograph (MPFS, [6]), with a spectral resolution AX ~ 8 A in 
4000-70001 spectral range. MPFS field consists of 16 x 16 1" x 1" spaxels. 
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For the nebula size (in its brightest part) 1" x l."5 it means rather poor 

spatial resolution, yet one can resolve the structure of the object, considering 
offset pixels, were the PSF wings of the closest parts of the nebula contribute. 

Data reduction system was written in IDL environment, using procedures, 
written by V. Afanasiev, A. Moiseev and P. Abolmasov. We also added at- 
mospheric dispersion correction in order to culculate correctly the nebula 
barycenters. 

Results 

More complete version of the results will be published in [7]. Integral fluxes, 
velocities, FWHMs and barycenters for selected lines are listed in Table 1. The 
integral unreddened flux in H/3 line is Fh0 = (2.00 ± 0.06) 10~"'^^erg cm^^s^^. 
Assuming Ha /Hp = 3 wc find interstellar absorption Ay = l.'"24. Unred- 
dening was produced using Cardelli et al. [5] algorythm. In the Table 1 we 
present the most important parameters of selected emission lines, including 
barycenter shifts versus Chandra source coordinates [2] . An interesting result 
can be seen that the barycenter is shifted to the West for low-excitation lines 
(Balmer lines, [SIT], [01]), and to the East for the most high-excitation line 
HeIIA4686. That can be understood as a global physical conditions gradient: 
the western part of MF16 is denser and colder than the eastern. 



Table 1. MF16 selected emission line properties. 



line 


F{\)/F{H0) F{\)/F{Hp) 
(unreddened) 


V, km/s 


FWHM,k Aac,arcsec 


A6c, arcsec 


HeII\A&m 


0.20±0.03 


0.22±0.03 


-5±18 


6.8±0.7 


0.302±0.05 


0.146±0.03 




l.OOiO.03 


1.00±0.03 


-28±4 


7.0±0.2 


-0.11±0.02 


-0.10±0.02 


[O/7/]A5007 


7.42±0.15 


6.94±0.14 


-17±2 


6.7±0.1 


0.058±0.015 


0.011±0.015 


[O7]A6300 


1.42±0.08 


0.93±0.05 


25±9 


10.1±0.4 


-0.196±0.015 


-0.237±0.015 


Ha 


4.73±0.09 


2.96±0.06 


-16±2 


7.8±0.1 


-0.111±0.013 


-0.08±0.013 


[Af77]A6583 


4.16±0.08 


2.59±0.05 


-19±2 


7.6±0.1 


-0.061±0.013 


-0.064±0.013 


[577]A6717 


2.46±0.02 


1.496±0.015 


-37.5±1.1 


7.7±0.1 


-0.116±0.014 


-0.081±0.014 



For some of the line profiles, that seem to have distinct two-component 
structure, wc also produce two-gaussian fit (see Table 2). The velocity shift 
between the line components is about 120/cm for all the lines resolved into 
two components. This value is consistent with the velocity gradient seen for 
[OIIIJA5007 line, exceeding 50knis~^ [7] and directed along the major axis 
of the nebula. Recent observations with higher spectral resolution [8] confirm 
the result, that the nebula expands with either ~ 120 (if we suggest just one 
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Table 2. MF16 two-component emission lines properties 



A,A Fi{\)/F{Hfj) Vi,krn/s FWHMi,k F2{\)/F{Hfj) V2,km/s FWHM2,k 

HellMQm 0.107±0.007 -142±20 11.6±0.6 0.093±0.004 10.±5. 4.26±0.15 
Hp 0.424±0.037 -122±17 10.24±0.50 0.58±0.02 -2±5 5.67±0.16 

[Om]A4959 0.875±0.014 -122±3 8.5±0.1 1.486±0.012 8±1 5.78±0.12 
[O7/7]A5007 2.696±0.015 -115±1 8.77±0.04 4.455±0.012 11±2 5.63±0.05 



shock wave with precursor), or ^ 60A:tos~^ (if we sec the both rececding 
and approaching parts of the nebula) velocity. In all the profiles the blue 
components are significantly broader, corresponding to velocity dispersion 
300-400fcms-^ 

Having the spectral line fluxes and kinematics, one can judge about the 
nature of the nebula. It is well-known [10, 2] that the optical emission-line 
luminosity of MF16 exceeds very much that of a SNR with comparable size, 
so one should suggest either a more powerful shock or an additional pho- 
toionizing source. If we consider a pure shock-wave excitation, the flux 
density will be connected to the shock parameters as was stated by Dopita 
& Sutherland [11]: 

FhI3 = FHf3,shock + ^precursor = 7.44 X 10"'' (loOftms-0 ^ (^^) + 

9-86 X 10-6 (Toofcr)'-'' X (^) ergcm-^s-^ 

Here ni and n2 are pre- and postshock number densisties, correspondingly. 
Basing on this formula and assuming the nebula a sphere with 13pc radius, 
and 712 — 500cm-^ [10] one can see the potential shock velocity varies from 
150km (for m = lOOcm"^, adiabatic shock) up to 350km (for m 
Ocm~^, fully radiative shock). 



High-excitation lines 

[OIII] and Hell lines in MF16 appear to be unusually bright for a shock- 
ionized nebula. Observed ratio of [OIII]A5007 / Hp ^ 7 requires a shock 
velocity about 400fcms-^ [12]. As for Hell / H^ ~ 0.2 value, it cannot be 
explained by existing shock -|-precursor models [12, 11] 

Hell ionization suggests a quite specific photoionizing source - extreme 
ultraviolet (EUV) with A < 228A. An estimate of the EUV luminosity can 
be made using ionizing quanta number estimate [13]: 

Euv CUB T .,„39 _i 

'^a<228;^:b^^177 XWA4686 -10 evg s 

A4DS(3 a^gjj ;^4ggg 
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Photoionization Modelling 

We have computed a grid of CLOUDY96.01 [14] photoionization models in 
order to fit MF16 spectrum avoiding shock waves. We have fixed X-ray spec- 
trum known from Chanrda observations [2] assuming all the plasma is located 
at lOpc from the central point source, and introduced a blackbody source with 
the temperature changing from 10^ to lO^K and integral flux densities from 
0.01 to 100 erg 001-"^ s-'^ . 




1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 -2.0 -1.5 -1 .0 -0.5 0.0 

loqioC[NII]?*6583 / Ha) logio{HelM686 / H/S) 



Fig. 1. [O///]A5007/H/3 versus [iV//]A6583/i/a (a) and //e//A4686 / Ha versus 
[0///]A4363 / [O///]A5007 (b) ionization diagrams with the integral spectrum of 
MF16 shown by a cross of error bars. Photoionization grid descibed in the text 
is also shown. Solid lines trace the constant temperature levels, dot-dashed lines 
correspond to constant black body flux densities. Only models with \gT{K) = 
4.80 - 5.30 and F = 0.1 - T.Oergcm"^ are shown. 

The best fit parameters arc \gT{K) = 5.15±0.05, = O.eiO.lerg cto^^ g-i^ 
that suggests quite a luminous ultraviolet source: Luv = (7.5 ± 0.5) x 
lO^^er^s"^. More completely the modelling results and methodics will be 
presented in [8]. The UV source is more than 100 times brighter then what 
can be predicted by extrapolating the best-fit model for X-ray data [2]. 

Conclusions 

MF16 has a nontrivial omission line spectrum, similar to those observed in 
NLRs of Seyfert galaxies. The nebula is highly asymmetric, with a dense cold 
western part, possibly connected with the radio-bright region [15]. Kinemati- 
cal properties suggest an expanding non-spherical shell. The sources of ioniza- 
tion acting in different parts of the nebula are probably different, strong ra- 
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diative shock at the western end and in the outer regions, and hard UV/EUV 
radiation source in the inncr/castcrn regions. 

Basing on the HeII4686 luminosity we should suggest a ~ lO^^erg 
EUV (A < 228 A) source responsible for the second ionization of He. CLOUDY 
simulations suggest a L(jv ~ W^^erg FUV luminosity needed to produce 
the actual [OIII]A5007 / H/3 ratio. 

The research has been supported by the RFBR grants 04-02-16349 and 
03-02-16341 and also by the joint RFBR/JSPS grant 05-02-19710. 
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